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We study the dynamical evolution of the so-called chiral magnetic effect in an 
electromagnetic conductor. To this end, we consider the coupled set of corresponding 
Maxwell and chiral anomaly equations, and we prove that these can be derived from 
chiral kinetic theory. After integrating the chiral anomaly equation over space in 
a closed volume, it leads to a quantum conservation law of the total helicity of the 
system. A change in the magnetic helicity density comes together with a modification 
of the chiral fermion density. We study in Fourier space the coupled set of anomalous 
equations and we obtain the dynamical evolution of the magnetic fields, magnetic 
helicity density, and chiral fermion imbalance. Depending on the initial conditions 
we observe how the helicity might be transferred from the fermions to the magnetic 
fields, or vice versa, and find that the rate of this transfer also depends on the scale 
of wavelengths of the gauge fields in consideration. We then focus our attention 
on the quark-gluon plasma phase and analyze the dynamical evolution of the chiral 
magnetic effect in a very simple toy model. We conclude that an existing chiral 
fermion imbalance in peripheral heavy ion collisions would affect the magnetic field 
dynamics and, consequently, the charge-dependent correlations measured in these 
experiments. 
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I. INTRODUCTION 

Quantum chiral anomalies have a serious impact on the macroscopic physics of hnite 
systems made by (quasi) massless fermions. One specihc example of this fact is the so- 
called chiral magnetic effect (CME), which refers to the generation of an electromagnetic 
(EM) current proportional to the magnetic held in systems where there exists an imbal¬ 
ance between the populations of right- and left-handed fermions [1]. Initially discovered in 
electroweak plasmas [2], it was pointed out that the CME could be realized in noncentral 
heavy-ion collisions (HICs), where large magnetic helds are generated, and the topological 
properties of QCD allow for the generation of the chiral fermion imbalance [31II]- The CME 
can also occur in condensed matter systems, such as Weyl and Dirac semimetals with chiral 
quasiparticles UU- Experimental observation of the CME in a Dirac semimetal has been 
recently reported through measurements of its magneto-transport properties [7]. 

The CME has been reproduced with a variety of methods like quantum held theory Eh 
HIHI, holographic approaches EHm and lattice studies [I21E3], where as an initial condition 
one assumes a constant and homogeneous chiral fermionic imbalance in the presence of a 
magnetic held. 

A formulation of transport theory for chiral fermions —that also describes ehects of 
the quantum chiral anomalies— has been developed in Refs. [TTlfTS] . For example, this 
can be achieved by considering the hrst h-correction to the classical dynamics of chiral 
fermions. Quantum corrections to the standard classical Boltzmann equation can be easily 
incorporated, and the resulting framework not only describes correctly the chiral anomaly 
equation for hnite temperature and density fermionic systems, but also other anomalous 
ehects, such as the CME. The basic requirement of the chiral kinetic theory (CKT) is the 
existence of well-dehned chiral quasiparticles. 

In this manuscript we study the dynamical evolution of the CME by deducing the form of 
the electromagnetic (EM) current from the CKT within the relaxation time approximation, 
for times larger than the collision time. Then, we study the Maxwell equations which have 
this EM current as a source, coupled to the chiral anomaly equation, which is also deduced 
from CKT. Similar coupled anomalous Maxwell equations (AMEs) have been hrst considered 
for small wavelengths in cosmological contexts; see e.g.. Refs. [T91 - I22] . Here, they are derived 
in a more general way. 
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The chiral anomaly equation expresses a quantum conservation law. After integrating it 
over space in a closed volume, it relates the helicity of fermions with the Abelian Chern- 
Simons number or magnetic helicity (as known in magnetohydrodynamics) in such a way 
that a linear combination of these two quantities is conserved. In the main part of this 
manuscript —regardless the method used to derive the AME— we study the dynamical 
evolution of the magnetic fields and of the magnetic helicity, constrained by this quantum 
conservation law imposed by the chiral anomaly. We assume that other processes that 
might change the fermion chirality are absent, or at least negligible in the range of time 
scales we consider. Depending on the initial conditions we observe how the helicity might 
be transferred from the fermions to the gauge fields, or vice versa, and hnd that the rate of 
this transfer also depends on the scale of wavelengths of the gauge fields in consideration. 

While the above features are generic for any chiral fermionic electromagnetic conductor, 
we focus our interest in the quark-gluon plasma (QGP) phase that is supposed to be realized 
in HICs. We employ a very simple toy model for the QGP, but that can be very illustrative 
to extract relevant conclusions. In particular, we assume a QGP in a static finite volume, 
and consider the evolution of an existing chiral fermion imbalance in the presence of a large 
initial magnetic held, with Gaussian shape, and no initial magnetic helicity. We observe 
a rapid creation of magnetic helicity and a reduction of the fermion chiral imbalance, of 
oscillatory behavior at very short time scales. We claim that the specihc evolution of both 
chiral imbalance and magnetic helicity depends on the initial spectrum of magnetic helds, 
but the effects we found are quite generic and have been overlooked so far in studying the 
impact of the GME in HIG physics. Further, we also consider the generation of a chiral 
fermion imbalance from magnetic held conhgurations with high magnetic helicity. 

This manuscript is structured as follows. The set of AMEs is studied in Sec. [TTj We rewrite 


the AMEs in Fourier space in Sec. II A , making use of circular polarization vectors, as then 
the magnetic helicity can be written in terms of the left- and right-handed polarization 
components of the magnetic held. We show how the dynamical evolution of each held 
component is diherent in the presence of a fermion chiral imbalance. Exact solutions for 


small wave numbers are found and displayed in Sec. IIB where we also show plots of the 


dynamical evolution of the axial chemical potential and the magnetic helicity. In Sec. IIG 


we treat the case of large wave numbers, where explicit solutions of the AMEs can only be 
numerically obtained. We point out that solutions of these equations show a very diherent 
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behavior for large or small wavelengths. In Sec. we focus on the QGP, and argue that 
even when the non-Abelian gauge held contribution to the chiral anomaly equation should 
be considered, at late times the color dynamics can be safely ignored. We present numerical 
results for the evolution of the chiral fermion imbalance and magnetic helicity assuming an 


initial Gaussian spectrum for the magnetic held. We present our conclusions in Sec. |IV[ In 
Appendix we give a quick review of the GKT and from it, derive the form of the EM 
current for a hnite temperature plasma with an initial fermion chiral imbalance, using the 
relaxation time approximation. In Appendix we give details on how we have numerically 
solved the AME. Finally, in Appendix we compute the helicity-hipping rate of Gompton 
scattering and argue that this process, together with the helicity change by Landau damping 
can be ignored in the context of HIGs. We use natural units {h = c = ks = 1), throughout 
the manuscript. 


II. ANOMALOUS MAXWELL EQUATIONS 


In this section we focus on the set of dynamical equations that are derived from the 
chiral transport approach. In essence, we study Maxwell equations with the EM current 
associated to the chiral fermions, together with the chiral anomaly equation, which is a 
direct consequence of the chiral transport equation. 

Similar equations have been used for electroweak plasmas in cosmological contexts in 
Refs, [ini EHl [22] by assuming the same form for the EM current as the one obtained here. 
For an electromagnetic conductor, it is possible to obtain a dynamical equation for the 
magnetic held. To this end we consider the Maxwell equation 

f + J = VxB. (1) 


where the EM current contains both Ohm’s law and the GME as derived in Appendix 


J 


Ns 

aE + J2 

S=1 




B . 


( 2 ) 


We take the curl of Eq. Q, and then use Faraday’s law and the explicit form of the current 
in Eq. (|^ to reach to 
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where we have defined the factor 






( 4 ) 


S=1 


that takes into the account the sum of the charges squared of every fermion species in the 
plasma and a = e^/(47r) is the EM fine structure constant. 

To arrive to Eq. ([^, we have assumed that both the axial chemical potential ps and the 
electrical conductivity a are homogeneous, or at least, that their spatial variations can be 
neglected in front of the spatial variation of the magnetic field. Otherwise, gradients of these 
two quantities would correct Eq. (|^. 

The magnetic field evolution is then linked to the value of the chemical potential of the 
chiral imbalance. Thus Eq. (|^ should be studied coupled to the chiral anomaly equation. 


Eq. (A6). In a system with different fermion species, we integrate the total axial current 


over the volume V to get 

d{nR — hl) 2Ca 1 


dt 7T V _ 

where Ti is the magnetic helicity density, defined as 


d^xE-B = 


Ca dTi 

TT dt 


(5) 


'H(t) = ^ [ d^xA 
k Jv 


B . 


( 6 ) 


This quantity is gauge invariant provided that B vanishes at, or is parallel to the boundary 
of V. As we will focus our study on closed systems, the use of the magnetic helicity as 
a dynamical variable will not bring any gauge dependence in our results. If we further 
take into account that there might be other processes that change the fermion chirality, the 
dynamical evolution of the chiral imbalance is then governed by the equation 

Cadn 


dn^ 

dt 


71 dt 


T/Ug , 


(7) 


where the helicity-fiipping rate T f depends on the system under consideration. This coeffi¬ 
cient will be extensively discussed later for a QGP. 

It is also convenient to express the above equation in terms of the axial chemical potential. 
At very high temperature T pg, one writes 

It = 

where ca = and yg = is the susceptibility of the chiral charge. 
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The interesting reason for writing the integrated form of the chiral anomaly eqnation in 
terms of the time derivative of the magnetic helicity is that it expresses a conservation law. 
In the absence of processes that change the fermion chirality, i.e. when Tj = 0, then Eq. Q 
expresses that the qnantity X = ns + is conserved. This means that an initial fermion 
chiral imbalance can be converted into magnetic helicity, but also that a magnetic held with 
helicity might be converted into a chiral fermion imbalance. The last mechanism has been 
recognized in cosmological contexts [19], but we point out that it is a generic feature of chiral 
systems that has been unnoticed in many studies of the CME, where a time-independent 
value of /is is generically considered. 

Note that if Tj = 0 the chiral anomaly equation allows one to make an estimate of 
the maximum magnetic helicity density that might be created in the presence of an initial 
fermion chiral imbalance as 

■Hmax ~ 

In practice, for the cases we will consider, this maximal value is not always achieved, as the 
magnetic field dynamics enters into a dissipative regime that does not allow for a sustained 
growth of the magnetic helicity. 

In the remaining part of the manuscript we will study in more detail the anomalous 
Maxwell equations and consider how the helicity might be transferred from the fermions to 
the gauge helds and vice versa at time scales f <C l/Tj. A different solution to the chiral 
anomaly equation has been considered at large times, in the regime t S> l/Tj [^, and in 
the presence of external parallel electric and magnetic helds. Such a solution allows one to 
hnd a peculiar magnetic held dependence of the electrical conductivity tensor that has been 
observed in Dirac semimetals [7|. In this work we will not consider the presence of hxed 
background external helds, and rather concentrate on the evolution of the chiral plasmas at 
shorter time scales. 


A. Solving the AMEs in Fourier space 

It is convenient to solve the AMEs. (|3]-[^ in Fourier space. After dehning the orthonormal 
set of vector polarization vectors (e+, e_, k) describing circular polarized waves [23] we dehne 
the Fourier mode associated to the vector gauge held potential 

Ak = -|- A^e_ + A^k . 


( 10 ) 
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If the vector potential is real then 

(4)' = , 

and the magnetic held in Fourier space is then expressed as 

Bk = -ik X Ak = -Bke+ + -Bk^- = -k{A^e+ - A,^e_) . 


( 11 ) 


( 12 ) 


It is possible to express the magnetic helicity density in terms of the circular polarized 
components of the magnetic held as 




(iflk 
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(13) 


where dflk is the solid angle element, not to be confused with the Berry curvature. 

Written in this form, we observe that the magnetic helicity gives account of an existing 
asymmetry in the magnetic held spectrum, measuring the diherence between its left- (-I-) 
and right-handed (—) polarization components. Note also that when 77 is expressed in this 
form one clearly observes that the chiral anomaly equation ([^ relates the imbalance between 
left- and right-handed helicity components of both fermions and magnetic helds in a quite 
parallel way. 

We note that the magnetic energy density can be written as 


pB — 
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Using circular components, the coupled equations ([^ and Q reduce to 


1 dB^ 


dn^ 

dt 


Ca 


a OB 

dflk 


+ 


dt 


a 


CaiJi^k\ 


B 


k ) 


dvr 






TTcr / 


(14) 

(16) 

(16) 


We see that the main ehect of the existence of a fermion chiral imbalance is that the left 
and right circular polarized components of the magnetic held follow a diherent dynamical 
evolution. It is easy to understand that in the presence of an axial chemical potential a 
nonvanishing magnetic helicity soon arises as a consequence of this wave asymmetry. 

In the remaining part of this section we will work in the limit ps T and use Eq. ([^ 
with a Ca not depending on ps. Hence, we can obtain simple solutions and estimates for 
the dynamical evolution of the chiral imbalance and magnetic helicity. However, in all our 
numerical examples we will consider the most general scenario, that may allow for larger 
values of the axial chemical potential. 










B. Solutions for small frequencies (or wave numbers) 


For small frequencies the coupled system of equations described in Eqs. (15 16) greatly 
simplifies, and exact analytical solutions can be found. In this case, the equations have 
already been discussed in the high-T limit (where ca oc does not depend on /is) to study 
the magnetic held dynamics in cosmological contexts 1201EU — when one is interested in the 
behavior of the system at large distances. These references focused on the appearance of a 
chiral magnetic instability, that might explain the existence of large cosmological magnetic 
helds. Here, we review the form of the solutions, and explain why a chiral magnetic insta¬ 
bility occurs. In addition, we also present other solutions without instability that show a 
growth of the magnetic helicity as a consequence of the chiral fermion imbalance. Further, 
we show that an initial large magnetic helicity might create a chiral fermion asymmetry. 

Let us assume that the system is isotropic. If we neglect the second-order time derivative 


in Eq. (15) by assuming frequencies much smaller than the EM conductivity a; -C a, then 


the solution of this equation reads 


exp 


a 


(17) 


where is the initial condition for the magnetic held, and we have dehned 


K, - ^ 




'0 


TT 


(18) 


We note that when /is is time independent then Kp = 


The value of /is is then found by plugging the solution for the magnetic held above 


in Eq. (16) and solving the resulting integro-diherential equation. The solution can be 

fJ-bit) = fi 5 ,o exp (-F ft) - CACa[H{t) - 7/(0)] , (19) 


expressed as 


and depends on the initial values of the axial chemical potential /is,o = h 5 (t = 0) and 
magnetic helicity 7/(0), and on the chirality-hipping rate F j. 

We will now review, in a general setting, how the chiral magnetic instability arises. This 
was hrst discussed in Refs. [201 EU, assuming an initial nonvanishing small value of the 
magnetic helicity, and in the approximation Fj 0. Let us hrst assume that the magnetic 
held conhgurations are maximal helical i.e. = 0. Then, the system of equations greatly 
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simplifies. Further, if one assumes a monochromatic magnetic helicity such that 




( 20 ) 


then the helicity at time t can be expressed as 

1-L{t) = ho exp 

where ho = 


— {Kp - ko) t 

a 


( 21 ) 


rf/ts 2c^Cako / Ca , 

-^1—As-i,|W(i). 


( 22 ) 


Then, one sees that for ko > Kp the magnetic helicity decays exponentially, as the Ohmic 
dissipation dominates its large time behavior. However, if ko < Kp the magnetic helicity 
grows exponentially, signaling the presence of an instability. The chiral imbalance then 
evolves until it reaches a tracking solution for p-s, 

Tlko 


h5,tr — 


Ca ’ 


(23) 


at which the chiral imbalance remains constant in time. In this situation the growth of 
the magnetic helicity is exactly compensated by its decay due to the Ohmic dissipation. 
Therefore, as K remains constant, ps does not change in time. 

The time scale where the instability starts to become effective can be estimated from the 


solution in Eq. (21) as 


t 


a 


'KCT 




2ko{Kp — ko) 2fcoC'a(p5^o ~ h5,ir) 


(24) 


and it clearly depends on how far the initial value of the chemical potential is from the 
tracking solution. 

Looking at the equation for ps, the characteristic time scale for the decrease of ps has a 
different value 


vrcr 


Tia 


2/coC'^a^CAa'H 2koC‘^a‘^CAho ’ 


(25) 


which also depends on the initial value of the magnetic helicity density. 

We show a particular example of the chiral magnetic instability in Fig. [T] In this re¬ 
alization we set an initial ps^o = 150 MeV and use maximal helical magnetic helds (with 
-®ko ~ !*)• make the arbitrary choice ~ where |Bo| is the magnitude 
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FIG. 1: Chiral magnetic instability. We use the parameters T = 225 MeV, a = 0.0244T, and 
consider a sphere of radius R = 10^ fm, ko = 0.5 MeV, /rs^o = 150 MeV, Fj = 0, and initial magnetic 
field as specified in the text. Left panel: Chiral imbalance as a function of time (black solid line) 


and tracking solution given by Eq. (23) (blue dotted line). Right panel: Magnetic helicity density 


as a function of time (black solid line), analytic solution at early times from Eq. (26) (red dashed 


line), and maximum value for the magnetic helicity density as given by Eq. (27) (blue dotted line). 


of the initial constant magnetic field in configuration space, that we take as |eBo| = 
where = 135 MeV is the pion mass. The other relevant parameters to determine the 
dynamical evolution are given in the caption of Fig. 

In the left panel of this hgure, we observe the decrease of the initial chiral imbalance 
according to the AMEs with ko < Kp. The dynamics is frozen when the chiral imbalance 


reaches the tracking solution in Eq. (23). In blue dotted line we show the value of this 
estimate. 

In the right panel of Fig. we present the magnetic helicity density (called simply 
“helicity” in the hgures) as a function of time. In black solid line we plot the numerical 
solution as obtained from our code. At initial times {t a/2kl) an analytical estimate of 
the magnetic helicity can be obtained from Eq. 

2koCa 


H = hn 


1 + 


an 


(h5,0 — h5,tr) t 


(26) 


This linear behavior is plotted using a red dashed line in Fig. When the tracking so¬ 
lution for the chiral imbalance is reached, the magnetic helicity density remains constant. 
The maximum value of the magnetic helicity density can be determined from the anomaly 
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equation Q with F/ = 0. It reads 

^max 


TT 


— ^0 + (^ 5,0 — n^,tr) 1 


(27) 


where and n^fi are the tracking and initial value of the chiral number density, respec¬ 
tively. The value for Tfmax is plotted in the dotted blue line in the right panel of Fig. 
providing an additional check of our numerical code. 

In the presence of several Fourier modes, the authors of Ref. [20] realized that there is 
an interesting inverse cascade phenomenon, where the magnetic helicity is transferred from 
the highest to the lowest modes, while the fermion chiral imbalance is washed out if one 
includes a small and nonvanishing value of F f. For a particular consideration of this case, 
we refer the reader to Ref. 1201 . 

We consider now a different situation, where the magnetic helicity might grow without the 
existence of an instability. As before, let us hrst consider a monochromatic initial magnetic 
helicity in Fourier space, such as 


Kf ^ - h) , 


±| 2 ; 


(28) 


but with |i?o f = |i?o'r = l-^on so that one has an initial vanishing magnetic helicity 
77(0) = 0. From the explicit solution of the dynamical evolution of the magnetic field, one 
can write the magnetic helicity density at time t as 


77(f) = 2hoexp 




a 


sinh 


2koKp 
-i 

a 


(29) 


where ho = ■ In the presence of a chiral imbalance, the magnetic helicity grows at 

short times, although at large times it decays exponentially, due to the Ohmic dissipation. 
More specihcally, for 7 <C ^ 

Ah rv 11 ^ n 

(30) 




na 

This case clearly differs from the one discussed previously, as in the presence of the chiral 
instability the growth is exponential, and it only stops when ps reaches the tracking solution. 
In this case, the maximal growth of 77 can be estimated as follows: the Ohmic dissipation 
takes effect at times 7 ~ By this time, the magnetic helicity will have reached the 
approximated value 
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2kl 


77 7 = —7T ~ ho sinh Ri ho sinh 


7F. 


ko 


V tt/co 


(31) 
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The growth of the magnetic helicity will depend on both the initial values of the magnetic 
held and chiral imbalance, but also on the value of k^. The growth can be very signihcant 
if fco < 

u ^ 

Quite similarly, one can obtain the time evolution of the magnetic energy density. One 
hnds 


Pb(Q = Pb{^) exp 


a 


cosh 


‘^kpKp 

a 


and its maximum growth —before the Ohmic dissipation takes effect- 


Pb{0) cosh ( ^ 


ps(0)cosh 


(32) 

can be estimated as 

(33) 


k„ J 7!ka 

Now, it is easy to realize that the increase of the magnetic helicity comes together with 


a lowering of the axial chemical potential, as dictated by Eq. (19). In the presence of a 
small nonvanishing helicity-hipping rate Tf, would also tend to vanish at later times. 
However if T/ 0, the chiral anomaly (which expresses the conservation of the helicity 
of fermions and gauge helds) dictates that the fermion chiral imbalance should grow when 
the magnetic helicity enters the regime where the Ohmic dissipation is effective. This is 
illustrated in Fig. [^for a plasma of hxed spherical volume, where we study both the fermion 
chiral imbalance and magnetic helicity density for a monochromatic field, in the regime 
k a. In all the examples in this section we assume ~ ^p^T'^, and make the arbitrary 


choice \Bn 


-, where |Bo| is the magnitude of the initial constant magnetic held 
in conhguration space, and we take |eBo| = where m-jr = 135 MeV is the pion mass. 
The explicit values of the other relevant parameters are detailed in the hgure’s caption and 


are taken from some of the conditions that are found in a QGP; see Sec. Ill 


In Fig. one observes a temporary creation of magnetic helicity due to the decrease of 
the initial chiral imbalance p^^ = 10 MeV. We also plot in the red dashed line the estimate 


Eq. (30) that approximates the increase of the magnetic helicity at initial times. 

Finally, let us describe the behavior of the system when several Fourier modes of the B 
held are taken into account. If we now assume 

N 


i=l 


ki) 


(34) 


where to ensure a vanishing initial magnetic helicity, we assume that iSo.q = l^oil — 
N is the number of Fourier modes. From the AME, all modes evolve independently 
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FIG. 2: Temporary generation of magnetic helicity at the expenses of an initial fermion chiral 
imbalance. We use the parameters T = 225 MeV, a = 0.0244T, and consider a sphere of radius 
R = 10^ fm, ko = 0.5 MeV, = 10 MeV, Fj = 0, and initial magnetic field as specified in the 
text. Left panel: Chiral imbalance as a function of time. Right panel: Magnetic helicity density 


as a function of time (black solid line) and analytic solution at early times given by Eq. (30) (red 
dashed line). 


and the total magnetic helicity density is the linear snm of partial helicities 

N N 


n{t) = = ^2 ho, iexp 


i=l 


2 = 1 


r 1 

—2—t 

sinh 


a 


a 


(36) 


47rQ: 


where ho^i = where we further assume ~ h 

Onr analysis shows that the dynamical evolution of the modes with higher ki would be 
much faster and would grow more than those of lower values of ki. At very short time 


scales, Eq. (30) tells us that the growth of Hi depends on the product hjho,*, and ho,i which 
depends on fcj. On the other hand, the helicity of the mode with higher k^ decays faster, as 
this happens at f ~ All the modes are coupled by the equation of the anomaly 


^ = A:.(^-Aw.(i). (36) 

Notice that when the first mode dies out due to Ohmic dissipation, the modes with a 
smaller wave number still contribute to the evolution of the chiral imbalance. Then, ps 
keeps evolving until the mode with smallest ki decays. This is illustrated in Fig. where we 
have used three Dirac modes ki = {0.2,0.4, 0.6} MeV plotted in red, blue and green lines, 
respectively, while the total magnetic helicity is plotted in black. Notice that according to 
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our model the amplitude of each mode ho,i depends on the wave number. The remaining 
parameters are chosen as in Fig. We show the estimate of Eq. (30) for each individual 
mode in dashed lines. 




time (fm) 


FIG. 3: Temporary generation of magnetic helicity at the expenses of an initial chiral imbalance 
for a polychromatic case. We use the parameters T = 225 MeV, a = 0.0244r, and a sphere of 
radius R = 10^ fm, ki =0.2 (red), 0.4 (blue), 0.6 (green) MeV, = 10 MeV, Fj = 0. Left panel: 
Chiral imbalance as a function of time. Right panel: In solid lines we plot the total magnetic 
helicity density (black curve) and individual contributions to it for each of the three modes (red, 
blue and green lines) as a function of time. In dashed lines we include the estimates of Eq. ( |30[ ) 
for early times. 


Finally, let us mention that Eqs. (15 16) also allow for a transfer of magnetic helicity 
density to the fermions. Interestingly, the chiral anomaly equation suggests that what¬ 
ever mechanism that may damp the magnetic helicity might be a source of chiral fermion 
imbalance. For an electromagnetic conductor, this mechanism is the Ohmic dissipation. 

A conhguration with an initial large magnetic helicity reduces its value, then creating 
chiral fermion imbalance. As an example, we present the solution for vanishing initial chiral 
imbalance but large magnetic helicity in Fig. For simplicity we use a single Fourier mode 
at ko = 0.5 MeV. The reduction of the initial magnetic helicity and the generation of chiral 
imbalance is maintained until the tracking value = T^ko/Ca ~ 107 MeV is reached. At 
this point both variables are constant in time, unless there is an independent process that 
reduces the chiral fermion imbalance, represented by a nonvanishing value of F/. This last 
case is represented in dashed lines, where we solve the same equation with a Fj = I/IO^ 
fm“^. 
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FIG. 4: Generation of chiral imbalance due to the presence of an initial magnetic helicity. We use 
the parameters T = 225 MeV, a = 0.0244r, R = 10^ fm, ko = 0.5 MeV, /rs^o = 0 MeV, Tf = 0. In 
dashed lines we plot the solutions to the AMEs when a Fj = 1/10^ fm“^ is included. Left panel: 
Ghiral imbalance as a function of time. Right panel: Magnetic helicity density as a function of 
time. 


C. Solutions for large frequencies (or wave numbers) 


While we have seen that analytical solutions can be found for the coupled system of 


Eqs. (15 16) when second-order time derivatives are neglected, in the most general case this 


is not possible, and a numerical study is required. An analytical solution is only possible if 
one further assumes that /is is constant. In such a case, assuming that the magnetic helds 
have a time dependence of the form ~ e*^^**, then the frequency is the solution of 


OJ-i 


1 2 p, 

iauj± — k ± -= 0 . 


TT 


For both left- and right-handed polarized components there are two solutions: 

icr . r—\ . . 4:k f, Ca^ 5 \ 

71 ) 


= y ± , 






(37) 


(38) 


dBi 


If we impose as boundary conditions B^{t = 0) = and 
for both polarizations read 


_ = 0, then the solutions 


= 2^k,o 


-V^) 


+ ( 1 + ^ 


(39) 


It is possible to check that when k a, Eq. (39) agrees with Eq. 10. for a constant /is 


This condition guarantees that oo a, which explains why the higher time derivatives can 
be neglected in the study of the system. 
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First, we observe that for large wave numbers a chiral magnetic instability can also occur. 


Let us consider that we deal with maximal helical helds, = 0. Equation (39) for can 


describe an exponentially growing solution whenever > 1- This occurs for values of 

the wave number satisfying k < , and assuming k > a. We also hnd that at the value 


of the tracking solution given by Eq. (23) the magnetic held remains constant. The time 
scale when the growth is effective is 

1 


U 


inst,h 


(\/^ “ 1 ) 


(40) 


a 


Note that for small wave numbers this occurs at a time scale given by Eq. (24). 


However, if we assume k ^ a and k > then both (j)± < 0, and the solutions for both 
left and right circular waves become oscillatory. 




cos 


aA±t 


1 . aA±t 

-— sm- 

A± 2 


(41) 


where we have dehned A± = \/\(j)±\. Then, we see that all the modes decay at time scales 
of order t ~ 2/(j, independent on the value of the wave number, although the oscillations 
at shorter time scales depend on its value. With this form of the solution, the magnetic 
helicity should decay at time scales of the order t 1/a. 

While we cannot present analytical forms of the solution if we assume a non constant /xs, 
the above solutions give us a hint of how the system behaves when we consider full time 
dependence of /xs. In this last case, the coupled equations have to be studied numerically. 
Details on how we have carried out such an analysis are given in Appendix 

We consider different situations with a fermion chiral imbalance for both monochromatic 
and polychromatic cases. In our numerical analysis we emulate the Dirac delta functions of 


Eq. (28) with a very narrow Gaussian, the amplitude of which is chosen to provide a correct 


limit to the Dirac delta function, see Eq. (Bl) 


We assume a Gaussian prohle for the magnetic held peaked at ko = 100 MeV and width 
K = 3 MeV, the amplitude is hxed as in the previous subsection. The magnetic helicity 
is chosen to be zero by taking the same value for the left- and right-handed polarized 
components of the magnetic held. The initial chiral imbalance is taken = 100 MeV. 
The results are shown in Fig. Being close to a Dirac delta with fco ^ and k > Ga/xs/vr, 
oscillatory solutions are found for and "H. These oscillations only disappear at t 3> 1/cr ~ 
40 fm. On average, the trend coincides with the hndings corresponding to the analytical 
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time (fm) 



time (fm) 


FIG. 5: Generation of magnetic helicity due to an initial chiral imbalance for a monochromatic 
field with ko = 100 MeV. As initial conditions, we set a chiral imbalance of /is^o = 100 MeV and 
zero magnetic helicity density. We use the same values of a and T that were chosen in previous 
figures, but reduce the volume of the sphere and take 7? = 10 fm. 


solutions for constant ^5 displayed in Eq. (41): the solutions are oscillatory and the initial 


chiral imbalance induces a temporary nonzero magnetic helicity until the latter decays at 
t ~ 1/cr. 

In Fig. I^we consider the same situation for a polychromatic magnetic held, with N = 3 
and wave numbers fcj = 100, 200, 300 MeV, and with the same specihc choice of dependence 
on ki of the initial spectrum that was used to analyze the evolution for small wave num¬ 
bers. There are interferences of the oscillations associated to every Fourier mode. Such 
interferences depend, however, on the initial wave number spectrum of the magnetic helds. 


> 

CD 



time (fm) 
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FIG. 6: Same plots as in Fig. but using a polychromatic field with wave numbers: ki = 
100, 200, 300 MeV. Notice the interferences among the three modes. 
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Finally, we consider the case in which a chiral imbalance is generated due to the presence 
of an initial magnetic helicity. We use the single Gaussian spectrum centered at = 100 
MeV and k = 3 MeV (monochromatic case). The initial amplitude of the magnetic held is 
taken as in the previous examples, but with a vanishing right-handed component. In Fig. 
we see how the initial magnetic helicity density generates a chiral imbalance which oscillates 
in time. Again, the helicity decays in times of the order of f ~ l/cr = 40 fm, creating a 
stable chiral imbalance that might eventually vanish by the effect of some F f. 



FIG. 7: Generation of chiral imbalance due to an initial magnetic helicity in the system. We set 
an initial vanishing chiral imbalance and a nonzero initial magnetic helicity density by choosing a 
left-handed-polarized magnetic field and a vanishing right-handed component. The profile of the 
initial magnetic field in Fourier space is a Gaussian centered at ko = 100 and width k = 3 MeV. 

Summarizing, we have analyzed the time evolution for monochromatic and polychromatic 
(and unrealistic) magnetic helds, and observe that the transfer of helicity from fermions to 
magnetic helds behaves very differently for large or small wave numbers. A much more 
realistic case should consider the full Fourier mode spectrum, and the hnal result would 
depend on the initial spectrum of magnetic helds. In the following section, we will show the 
behavior of the helicity transfer when the initial spectrum is a wide Gaussian for the QGP 
conditions. 


III. APPLICATION TO THE QGP 

It has been claimed that the GME might have several phenomenological consequences in 
relativistic heavy-ion collisions [siia^. In these experiments a chiral imbalance between 
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left- and right-handed quarks might arise as a result of an event-by-event fluctuation of 
the non-Abelian topological charge. While averaged over all events, the topological charge 
should be zero, a gluon held conhguration with nonvanishing winding number should be 
converted, thanks to the non-Abelian chiral anomaly, in a difference in the populations of 
left- and right-handed quarks. Further, in peripheral heavy-ion collisions, large values of 


the magnetic helds are expected to occur, of the order of e|Bo| ~ (1 — 10)m^ [25]. It has 
been claimed that the CME should have an imprint on the charge asymmetry huctuations 
detected in HICs [21|. Unfortunately, it is not easy to identify the occurrence of this quantum 
effect in these experiments, and different experimental signatures of the CME in HICs have 
to be studied. 

Studying the dynamical evolution of the chiral magnetic effect in the QGP is much more 
involved than in the systems discussed in the previous part of this manuscript. The main 
reason is due to the fact the chiral anomaly is also affected by the non-Abelian gauge helds. 


Thus, rather than Eq. (A6) the total axial current obeys 


= 


2Ca 


E B + 


asNf 


TT 


TT 


Ea B„ , 


(42) 


where and B^ are the chromoelectric and chromomagnetic helds, respectively, and is 
the QCD hne structure constant. We have neglected above (and will do it in the following) 
corrections to the chiral anomaly equation due to nonvanishing values of the quark masses 
rriq, which might be entirely justihed for high values of the temperature, if we are in the 
regime where rUg T. 

It is then clear that in the QGP the evolution of the quark chiral imbalance has to be 
studied together with both the Maxwell and Yang-Mills equations associated to the system. 
We should note that it has been claimed that in the presence of a quark chiral imbalance, 
there should be also chiral chromomagnetic instabilities, very similar to what happens in the 
Abelian sector [271 EH]. However, our claim is that this can only be so in the limit of very 
weak couplings, -C 1, and it is reasonable to expect that at large and also at moderate 
small values of ag, as those that are expected to occur in HICs, the nonlinearities of the 
gluon dynamics would saturate the growth of the chromomagnetic helds. This is at least 
what happens in the case of Weibel instabilities in a non-Abelian plasma [21]. This issue 
requires a further detailed study, that we will not address here. 

At the initial stages of the HIC, a nonvanishing value of the topological charge associated 
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to the gluon fields might create the quark chiral imbalance. At the very early stages of 
the collision, the chromoelectromagnetic fields are so strong that is reasonable to ignore the 
effect of the electromagnetic fields in the evolution of the quark chiral imbalance. In this 
situation the evolution of quark chirality might be affected by the contribution of sphaleron- 
like transitions. QCD sphalerons allow to change the Chern-Simons number associated to 
the gluon fields, and this change, due to the chiral anomaly associated to the non-Abelian 
gauge fields, leads to a variation of the quark chiral fermion number |30]. In particular, in 
Ref. [21], and making use of a fluctuation-dissipation result, the time evolution of the chiral 
quark number due only to QCD sphaleron transitions is expressed as 


dns _ _ {2Nff Tsph 

dt ~ X5 2T 


(43) 


where Fsph is the so-called sphaleron rate. 

The sphaleron rate has been computed in Ref. ED at very large T for a SU{Nc) theory, 
assuming not only a small value of the QCD coupling constant, but that the inverse of its 
logarithm should be small. This means that the results of Ref. [3T] are only strictly valid 
at very large values of T. The authors of Ref. [3T] suggested that if one could extrapolate 
their results to the regime where ~ 1/3 —a value more in accordance with what one 
expects to find in HICs— one could take Fsph ~ It is worth noticing here that 

the sphaleron rate has also been computed in a strong coupling regime using AdS/CFT 
techniques in Ref. [32], finding Fgph = Nc)'^T^/. If we assume that the sphaleron 
rate is correctly described by Fgph ~ assuming as = 1/3 one then concludes that 

the QCD sphalerons start to have an effect at time scales tgph ~ 0.45/T. Thus, even in the 
case that a chiral imbalance is created in the very initial stages of the HIC, as suggested in 
Ref. [23], one may expect that this evolves in time in the regime where color dynamics is 
relevant. 

The color charges and currents generated in the early moments of the HIC are rapidly 
bleached and have a shorter lifetime than the electromagnetic currents, so we expect they 
might be absent in the long time behavior of the system. In particular, in a weak coupling 
scenario the color currents are known to disappear at time scales of order t - ■ ^ 


Os log l/asT 

while for moderate to strong coupling one might expect bleaching to occur in shorter time 
scales. Thus, at longer times one can study the evolution of the quark chirality ignoring 
the effects of the the color gauge fields, and taking into account only the electromagnetic 
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fields. We will assume in the following that a chiral imbalance persists in the regime where 
there are no in-medium chromomagnetic helds. The axial charge created has been estimated 
in models where color flux tubes are generated in the early moments of HIC [351135] , and 
in the regime where color dynamics can be ignored, values of ps ~ 100 MeV can still be 
considered | 37 |. 

We analyze here the evolution of the chiral magnetic effect in a very simple toy model 
that would correspond to the time scales where the color dynamics might be ignored. Our 
primary goal is to evaluate whether the generic effects we have discussed so far —which arise 
in any electromagnetic conductor where its charged chiral fermions with chiral imbalance 
are coupled to the electromagnetic helds— might have any relevance in HICs, rather than 
presenting a much more realistic evolution of the chiral imbalance in HICs. Let us mention 
here that only Ref. [38] has considered the time evolution of the electromagnetic helds in 
the QGP in the presence of a chiral imbalance, however disregarding the chiral anomaly 
equation in studying the evolution of the magnetic helds. 

A. Toy model and numerical results 

In this subsection we study the evolution of the quark chiral imbalance and CME in 
a very simple toy model. We consider the presence of very strong magnetic helds, of the 
order of those found in noncentral HICs, with e|B| ~ [25|. We will ignore the ehects 

of the gluon helds as these have a shorter lifetime than those of the electromagnetic helds 
and currents. In particular, this means that we will ignore the sphaleron rate that allows 
us to change the quark chirality. We consider a hxed spherical volume V of radius i? = 10 
fm, made of an isotropic QGP. This choice of geometry might look too simple to describe 
the hreball resulting in a peripheral collision giving rise to large magnetic helds in the HIG, 
while a cylinder might look more appropriate. However, our analysis shows that the ehect 
we study does not depend on the geometry of the system. Although our toy model is of hxed 
volume and assumed to be inhnitely long lived, our only aim is to see whether the transfer 
of helicity from the fermions to the gauge helds can take place in time scales much shorter 
than the lifetime of the HIG hreball, as otherwise, it would be irrelevant for the realistic 
HIG. We also aim the order of magnitude of the ehect. We should warn the reader that a 
much more careful analysis valid for HIGs should consider the volume expansion and a more 
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precise spectrum of the magnetic fields which are obtained in noncentral HICs, which can 
be obtained with a study similar to the one carried out in Ref. |39j . 

For the QGP the chiral density is expressed as 


"5 = (l4 + , 


(44) 


where Nc = 3 and Nf = 3 stand for the number of colors and the number of light quark 
flavors, respectively. For simplicity, we assume that = (fr + hL)/2 T in all of our 
computations. Further, for three light quark flavors, C = iVc(| + | + |) = 2. 

We assume values of the chiral chemical potential of the order ps ~ 100 MeV, as suggested 
in Ref. [37], although different from that reference, we assume that it is homogeneous over 
the whole volume. We take values of the electromagnetic conductivity for QCD with three 
light quark flavors as found in the lattice QCD a = 0.0244 T [30], and we will assume 
temperatures as those found in the Relativistic Heavy-Ion Collider (RHIC) experiments, 
T = 225 MeV. This means that a ~ 5.4 MeV. 

We assume a Caussian prohle in Fourier space for the initial magnetic held, similar to the 
Caussian conhgurations taken in Refs. [371 E] to do studies of anomalous hydrodynamics 


for the HIC. The precise form of the magnetic spectrum is taken as in Eq. (Bl), similar to 
the one taken in the previous section as a regulator of the Dirac delta: 


Bkfi = bo exp 


1 f k — ko 

2 V K 


(45) 


In our toy model the peak of the Caussian is taken to be in the range between 100-300 
MeV, with a width of the order 40 MeV (consistent with the widths in conhguration space of 
Refs. [371131]). This means that if Fourier transformed to conhguration space, the magnetic 
held is concentrated in a small hnite region contained in the volume V. With this form 
of the initial spectrum we also guarantee that the typical wave numbers associated to the 
magnetic held are bigger than the inverse of the length scale of the size of the system ~ 1/R. 
As in this case we are not interested in taking the inhnitely narrow Caussian, we simply 
take the amplitude of the magnetic held in such a way that in conhguration space it takes 
the value |eBo| = Note that the allowed wave numbers of any magnetic held inside the 
volume have to be /c > 1/i? ~ 20 MeV, much larger than the value of the electromagnetic 
conductivity, fc 3> cr. 
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We solve numerically Eqs. (15) and (16) as described in Appendix [Bj Let us stress that 
several scattering processes between quarks and photons, quarks and gluons, and quarks and 
quarks give a contribution to the helicity-flipping rate Tf. However they scale as a^m?^/T 
and respectively, and it can be shown that they are relevant at time scales much 

larger than the expected lifetime of a HIC hreball ~ 10 fm (see Appendix for a much 
more detailed discussion). Because the relevant time scales associated to the change of the 
quark chirality due to the chiral anomaly are much shorter, it is reasonable to approximate 


r/~o. 

Let us hrst mention that the occurrence of a chiral magnetic instability is very unlikely, as 
this requires values of the magnetic held wavelengths k < . If we take the minimal value 

A; ~ 20 MeV, then the initial value of the chiral magnetic potential should be > 4304 
MeV, an extremely large value very unlikely to occur. For this reason, we ignore this scenario. 

However, a temporary creation of magnetic helicity is indeed possible in a HIC context 
if an initial chiral imbalance is present in the system. We assume a Gaussian shape for the 
initial magnetic held, centered at fco = 100 MeV and width of k = 40 MeV. We plot the 
results of both the axial chemical potential and magnetic helicity in Fig. |^up to large times, 
although relevant time scales for the HIC should be t < 10 fm. 



FIG. 8: Generation of magnetic helicity due to an initial chiral imbalance. The parameters are 
|eBo| = p 5 ^o = 100 MeV and a Gaussian initial magnetic field centered at ko = 100 MeV and 
width K = 40 MeV. 

We observe the oscillatory behavior of the chiral chemical potential with very small 
amplitude, of the order of 0.05 % with respect to the initial value, and the appearance of 
also an oscillatory magnetic helicity. These oscillations are characteristic of the solutions for 
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large wave numbers discussed in Sec. |II C[ although here we have a continuous spectrum of 
modes. 

Just to check the dependence of the results on the value of the magnetic field, we assume 
now the same situation with a value of the magnetic field as |eBo| = lOm^ but using the 
same parameters of the Gaussian as before and with the same chiral imbalance of = 100 
MeV. We show our results in Fig. The qualitative shape of the curves does not change. In 
particular, the periodicity of the oscillation remains invariant because it is solely controlled 
by the typical wave number, not by the strength of the magnetic field. In this case, the 
axial chemical potential presents larger fluctuations due to the bigger value of the magnetic 
field. We find variations on /xs around 4 % with respect to the initial chiral imbalance, that 
is, the fluctuations are 100 times larger than the previous example. The magnetic helicity 
density is also amplified by a factor of 100 with respect to the previous example (due to the 
fact that the magnetic field is chosen to be ten times bigger than before). 




time (fm) 


FIG. 9: Generation of magnetic helicity due to an initial chiral imbalance. The parameters are 
|eBo| = lOm^, = 100 MeV and a Gaussian initial magnetic field centered at ko = 100 MeV 
and width «; = 40 MeV. 


We now test the dependence of our results to the value of /cq. With initial ps.o = 100 
MeV and magnetic field |eBo| = we take = 200 MeV. The results are shown in 
Fig. 10 One can see an increase of the oscillation frequency with respect to Fig. This is 
because when increasing the typical wave number ko, the factor A± grows accordingly, and 
the sinusoidal oscillations have smaller periodicity. 

To conclude the study of our anomalous Maxwell equations we present the case with a 
dynamical generation of chiral imbalance due to an initial magnetic helicity in the system. 
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FIG. 10: Generation of magnetic Felicity due to an initial chiral imbalance. The parameters 
are |eBo I = = 100 MeV and an initial magnetic field taken to be Gaussian centered at 

ko = 200 MeV and width k = 40 MeV. 


In this case we impose = 0 MeV and a left-handed polarized magnetic field consistent 
with |eBo| = being the right-handed polarized component kept to zero in order to force 


a nonzero initial magnetic helicity. Our results are presented in Fig. 11 



FIG. 11: Generation of chirality imbalance due to an initial magnetic helicity. The parameters 
are: |eBo I = p 5 _o = 0 MeV and an initial magnetic field taken to be Gaussian centered at 
ko = 100 MeV and width k = 40 MeV. 


We observe a generation of chiral imbalance at the cost of the decrease of the initial 
magnetic helicity. The effect continues until times as large as 100 fm (or possibly earlier, 
depending on the presence of some nonzero Tf). Finally we study the same generation of 
chiral imbalance when the Gaussian prohle is centered at a higher wave number, ko = 200 
MeV. We show the numerical results in Fig. 
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FIG. 12: Generation of chirality imbalance due to an initial magnetic helicity. The parameters are 
|eBo| = = 0 MeV and an initial magnetic field taken to be Gaussian centered at ko = 200 

MeV and width k = 40 MeV. 

From the past examples, we expect an increase of the freqnency of the oscillations. In 
addition, we observe an increase of the generated chiral imbalance. This amplihcation is 
affected from the increase (a factor of 2) of the initial magnetic helicity. To nnderstand 


this factor we can compnte the initial magnetic helicity density snbstitnting Eq. (45) for 


the left-polarized magnetic held into Eq. (13) (and setting the right-polarized component to 


zero). Performing the integration and nsing Eq. (B4) we obtain 


«o - rr, —ko 


(46) 


an^V 

where the relation is valid for ko moderately larger than k. Therefore, the initial magnetic 
helicity is ronghly proportional to ko, which acconnts for the amplihcation factor of the 


magnetic helicity density in the last example. As opposed to the resnlts in Sec. IIB (where 
we neglected the second time derivative in the AME), the tracking chiral imbalance does 
not depend explicitly on ko, only indirectly throngh the initial magnetic helicity density. 
Therefore, we want to remark that a more precise knowledge of the initial magnetic held 
prohle is needed to determine a solid conclnsion of this ehect in HIC experiments. Note 
also that we have only shown the time evolntion of the magnetic helicity, and not of the 
single magnetic held components, which might evolve diherently depending on its initial 
geometrical conhgnration. A mnch more realistic conhguration for the HIC wonld have 
implied to choose the initial magnetic held created by moving charges in one direction, as 
done, for example, in Ref. |3H]. However, as the amonnt of helicity is governed by the valne 
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of /is and the chiral anomaly equation, we expect that our plots for "H would not change 
drastically. 


IV. CONCLUSIONS 

In this work we have studied the coupled system of Maxwell and chiral anomaly equations 
valid for EM conductors made up by chiral fermions. We have proven that a possible 
derivation of this coupled set of equations can be performed in the framework of chiral kinetic 
theory in the relaxation time approximation. The resulting set of anomalous equations are 
fully general, and can have a wide range of applications in systems so different such as 
electroweak and quark-gluon plasmas or in condensed matter systems such as Dirac and 
Weyl semimetals. 

We have considered the AME to study the evolution of both magnetic held and fermion 
chiral imbalance. After integrating the chiral anomaly equation over space on a hnite volume, 
a conservation law relating the chiral charge and magnetic helicity densities is obtained. 
When expressed in Fourier space and in circular polarized components, we observe that the 
dynamical evolution of the left- and right-handed components of the magnetic held in the 
presence of a hnite /is is diherent. Then, it is easy to understand that in the presence of 
a fermionic chiral imbalance a nonvanishing value of the magnetic helicity, which measures 
the asymmetry between the two diherent polarizations, soon arises. We have carried out 
a detailed study of the transfer of helicity from fermions to magnetic helds, with diherent 
initial conditions, and realized that the rate of this transfer depends drastically on the 
characteristic wave number of the magnetic held in consideration, while the amount depends 
on the magnitude of the initial magnetic held. 

Using the same set of equations but with diherent initial conditions, we also showed that 
a chiral equilibrated plasma might develop an imbalance between right- and left-handed 
fermion populations if a nonzero magnetic helicity is applied to the system. This occurs 
whenever there is a mechanism that damps the magnetic helicity, in the case we are studying, 
the Ohmic dissipation. Then the chiral anomaly equation dictates that the lowering of the 
magnetic helicity comes together with an increase in the chiral fermion imbalance. The 
equations tell us that the creation of fermion imbalance persists in the system as long as 
no other dynamical process erases it. We examined this case for systems with characteristic 



small and large frequencies, and found the rate of the transfers in each case. 

We have focused our attention to the QGP that is supposed to occur in heavy-ion col¬ 
lisions. While it has been recognized that the CME would have a clear impact on the 
charge-dependent hadron azimuthal correlations detected both in the RHIC and in the LHC 
for noncentral collisions, our main claim is that the CME changes the dynamical evolution 
of the magnetic held in a way that might affect the studies of those correlations performed 
to date A similar conclusion was reached in Ref. |38], after studying the magnetic 

held originated by a point particle in the QCP in the presence of the CME current, where 
also rapid oscillations of the magnetic held were obtained, even if the time evolution of ps 
was ignored. 

While we realize that the occurrence of a chiral magnetic instability is very unlikely 
to occur in HICs izg — as it would require extremely large values of the chiral chemical 
potential— a hnite nonvanishing magnetic helicity would soon arise if a fermion chiral im¬ 
balance persists in the late stage evolution of the QCP. The magnetic helicity gives account 
of a rather nontrivial topology of the magnetic held lines, with formation of knots and heli¬ 
cal structures. While our conclusions have been extracted from a simple toy model, a much 
more careful analysis should be carried out for applications to HICs. Then, another relevant 
consequence of the the chiral fermion imbalance is that it would ahect the estimate of the 
lifetime of the magnetic helds in a HIC. This is so because with magnetic helicity a magnetic 
held conhguration tends to be more long lived |12] . 

The possibility of a chiral imbalance generation due to an initial presence of helical helds 
also deserves more attention in the context of HICs. This case requires a deeper knowledge 
of the topological arrangement of the initial magnetic helds created by spectator nucleons. 
Any diherence on the right- and left-handed polarized components of the magnetic held 
might potentially cause nontrivial structures as knots, that induce the creation of a quark 
chiral asymmetry. It is important to recognize that this generation of chiral imbalance has 
an Abelian origin and could represent an alternative scenario to the one proposed so far, 
which is intimately linked to the topological properties of a non-Abelian gauge theory such 
as QCD. 

We leave for future work a more realistic study of the anomalous magnetohydrodynamic 
equations in the context of HICs. For this goal, one can follow the lines of Refs. 
where the particles of the HIC are described by a dynamical quasiparticle model. The 
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electromagnetic fields are computed from the Lienard-Wiechert potentials for all charged 
quasiparticles (both spectators and participants). Likewise, these retarded EM helds pro¬ 
duce a backreaction on the quasiparticles which are propagated in time. In this way, the 
spacetime evolution of electromagnetic helds can be computed with a realistic expansion of 
the hreball. The extension of this model to include a chiral imbalance will be addressed in 
the future. 
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Appendix A: Chiral kinetic theory in the relaxation time approximation 


In this section we show how the AMEs studied in the remaining part of the manuscript 
can be derived from CKT. First, we will briehy review the basic inputs of the (semi)classical 
kinetic theory (see Refs. [HHIHIIII! for a more detailed description). 

The distribution function associated to a chiral fermion fp obeys the transport equation 


dfp dfp dfp dfp 


(Al) 


where C[fp\ is the collision term. The fermion equations of motion take into account the 
effects of the Berry curvature and for a fermion of charge e read 

(1 + eB ■ r2p)f = V -|- e(E x Qp) + eB(v • fip), (A2) 

(1 -I- eB ■ f2p)p = e(E -|- v x B) -|- e^f2p(E ■ B), (A3) 

where we have dehned E = E — - -tA and v = ^ and the Berry curvature is expressed as 

e or op ^ 


nr. = H- 


p 


2|p| 


(A4) 


where H = ±1 gives account of the helicity of the right- and left-handed fermions, respec¬ 
tively. The energy for particles/antiparticles is also dehned as [TOl I 


Cp = ±p(l - eB ■ 


P/ 5 


(A5) 
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respectively. 

Then it is possible to show that the axial current defined as obeys the chiral 

anomaly equation [HHISIIIZ] 

■ B , (A6) 

while the vectorial current jy = is conserved. 

The system of equations presented here has been analyzed in the collisionless limit using 
a linear response analysis in Refs. [I61[I7]. There, it has been checked that the resulting EM 
current reproduces both the nonanomalous and anomalous hard-thermal-loop/hard-dense- 
loop effective actions that appear in a quantum field theory approach. 

The next natural step is to see the effects of collisions in the dynamics. So far, collision 
kernels in chiral kinetic theory have not been constructed, although it has been pointed out 
that they should be nonlocal in order to respect Lorentz invariance |18]. However, in order 
to see the effects of collisions at large times, one typically uses the so-called relaxation time 
approximation (RTA). Although the RTA gives a good qualitative behavior of the solutions 
of the Boltzmann equation at large times, it is well known that it does not respect the 
conservation laws of the system. In Ref. |19] an improvement to the RTA collision term 
was proposed, such that it does not spoil particle number conservation. We propose this 
alternative form for the relaxation-time collision term in CKT, not to violate vector current 
conservation, nor the chiral anomaly equation. Hence, we will use a collision term of the 
Bhatnagar-Gross-Krook (BGK) type [IH] 




1 

T 



n 



req 

Jp 


(A7) 


It is easy to prove that / + eB ■ flp)C^^^[fp] = 0, and thus that particle number 

conservation is not affected by the presence of the collision term. We note here that in 
Ref. [SU] a different form of the collision term was used. The authors checked a posteriori 
that the solution they found was respectful with the particle conservation laws. We note 
that the proposal we present above gives the same hnal answers, but in a much faster way. 

To study the effects of the collisions in the dynamics we proceed as in Refs. jlGlIIT]. The 
distribution function is expanded in powers of e up to linear order 


U = /” + e (■’/h + Vr’) 


(A8) 
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and the transport eqnation is solved in a pertnrbative expansion, assnming that the vector 
gange helds are of order = 0{e), while the spatial derivatives are d/dx^ = 0{6). Inserting 
this expansion into the kinetic eqnation, the relaxation time modihes the expression of the 
ont-of-eqnilibrinm distribntion fnnctions. The solntion can be easily compnted in Fonrier 
space 




E ■ V 


C'' 


.-1 


) = H, 


iv ■ k + dp 

V ■ 1 dC 
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iv ■ k + 
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n. 
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eq 


(A9) 

(AlO) 


iv ■ k + T ^2p dp 

where is the particle density compnted np to order 0(e). This solntion takes the same 
form as expressed in the collisionless limit (see Eqs. (36) and (38) of Ref. [T7]) bnt with 
the replacement iv ■ k ^ iv ■ k + in the denominators of the expressions, pins one new 
additional term 


&=-Un“‘ 

T 


n 




dklr 


dvr iv ■ k + T 


-1 ’ 


(All) 


where dQ^ is the solid angle element. 

Ont of eqnilibrinm, we parametrize the right-/left-handed particle densities in a similar 
way as in thermal eqnilibrinm. 


^ i^R/L , (A12) 

where we have introdnced ont-of-eqnilibrinm chemical potentials. Dehning = p — we 
then express — Sp = yd/i. It is then possible to show that the solntions 

we have fonnd here agree exactly with those of Ref. [50]. 

In a self-consistent treatment of the system, the EM helds which appear in the transport 
eqnation are generated by the same charged fermions. Thns, they shonld obey Maxwell 
eqnations 

= r . (A13) 

Assnming that the system is formed by np to Ng different particle species, the total electro¬ 
magnetic cnrrent is expressed as 

Ns 

= (P) J) = [(n,,j5) - (hsjs)] , (A14) 

S = 1 

where is the charge associated to a given species of particles, and for a relativistic plasma 
we have inclnded both fermions and antifermions with different chiralities of different species. 
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We are interested in the solutions in the limit where iv ■ k t Then one finds the 
total EM current 

Ns 2 

= + (A15) 

where jj,^ = Hr — HL) = rm'jj/S is the electrical conductivity of the plasma, = 

Debye mass, and we have summed the contribution of every 
fermion species. Thus, we see that the EM current at large times gives account of both 
Ohm’s law and the CME. 


Appendix B: Numerical study of the AMEs 


We present here some details on how we have performed the numerical study of the AME 
equations. Let us recall that for k ^ a one needs to consider the second time derivative in 


Eq. (15). However, in this case an analytic solution analogous to Eq. (17) does not exist for 


a nonconstant ps. The system of equations (15 16), together with a relation between and 


Ps, such as the relation (44) for the QGP, must be solved numerically to obtain (t) 


and We discretize the time, U^U + At, and apply the following scheme to solve 

the system: 

1. We set initial conditions for the chiral imbalance and the magnetic fields B^{ti). 

The latter is taken to have a Gaussian spectrum peaked at ko 


B^{ti) = 6oexp 


1 k k — ko 

2 V K 


(Bl) 


where the width of the Gaussian is controlled by k. 

In this work we have used two different criteria to fix the value of bo'- 


In Sec. IIB we solve the AMEs in the low wave number limit using a Dirac delta 
as an initial condition for the magnetic field squared. To be consistent with this 
choice in the high wave number case we use 


bo = 



B„|l/ . 


(B2) 
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This value allows us to connect the Gaussian profile with a Dirac delta spectrum 
in the limit 

2 


lim = fcolBol V" lim 


exp 


k — kr 


K 


= - ko) , 

(B3) 


in consistency with Eq. (20) of Sec. IIB 


For the QGP we do not attempt to take the Dirac delta limit. Therefore, we use 
more realistic (wider) Gaussians instead. For this reason, we just fix bo in such 
a way that performing an inverse Fourier transform we recover the value |Bo| in 
configuration space. Therefore, 

|Bo|(27r)3/2 


bo = 




(B4) 


for all cases in Sec. IIII Al 


2. We also need to fix the initial condition dB^ / dt{ti). To avoid an initial impulse of the 
magnetic helicity we set it to zero for all wave numbers. One can compute the value 
of the magnetic fields at ti + At as 


B^{ti + At) ~ B^{ti) + At ■ 


OB, 


± 


dt 


iU) 


(B5) 


3. Then, one transforms Eq. (15) into a finite-difference equation using the approximation 

B^{ti + 2At) - 2B^{ti + At) + B^{ti) 


d^B^ 

- 

This allows us to obtain the magnetic field at ti + 2At for all k 

,2 _ CaiJ,5{ti)k 


(B6) 


B^{ti+2At) ~ (2 — crAt) B^{ti+At) 


1 — a At + { k^^ ^ 


TT 


(At)^ 


This value is then used to extract the derivative of the fields at ti + At. 

^ ^ + 2Ai) - BiitO 


dt 


2At 


Bk (*i) ■ 
(B7) 


(B8) 


4. At ti we compute the rhs of Eq. (16) by numerical integration of the magnetic fields 


and their derivatives. With this, we obtain dnr,/dt{ti) to predict 

+ At) n5(i,;) + At ■ 


(B9) 
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5. Using Eq. (44) one can convert n^{ti + At) into + At). 


6. We iterate the last three steps until the hnal time t/. 

We need to probe several orders of magnitude in time. For this reason, we use a dynamical 
discretization step At that increases with time to speed up the calculation. 

Appendix C: Helicity flip scattering 


In this Appendix we justify the reason why one can neglect the helicity-flipping rate E f 
that enters into Eq. ([^ for the QGP in the time scales we are considering the system. 

There are several scattering processes that allow us to change the quark helicity. Let us 
first discuss Compton scattering in detail, when a massive fermion might change its helicity 
due to Compton scattering ESI- For almost massless fermions, helicity becomes synonymous 
with chirality. For simplicity, let us consider the situation where fiR, fiL -C T, and compute 
the rate in QED. We compute the helicity-flipping rate associated to Compton scattering as 

1 dn^ 6 / duR drii 


Tf = 


ns dt 


\ dt 


dt 


(Cl) 


neglecting completely the effect of the quantum anomaly. We compute the rhs of the Eq. (Cl) 
with the use of the Boltzmann-Uehling-Uhlenbeck equation for right-handed fermions 

where F denotes the combination of distribution functions 
F = /fd(P )/BE(fcO[l + fBE{k)][l - / fd ( p )] - /BE(fc)/FD(p)[l - /fd(P )][l + fBE{k')] (C3) 

and /fd and /be stand for the Fermi-Dirac and Bose-Einstein thermal equilibrium distribu¬ 
tions corresponding to the electrons and photons, respectively We integrate over p to obtain 
the rate duR/dt, 

rln rt r rl'^rt 

(C4) 


R 


d?p dfR 


dt J (27r)3 dt 

In the computation we ignore all the fermion masses in the kinematic factors, e.g. E{p) ~ p 
and concentrate on the leading term, in the scattering amplitude squared, which reads 


m^s 


\M - - I = 4e 

I-" ^ (_^ + ^2)2 ’ 


(C5) 
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where the mass in the denominator is kept to regularize the infrared limit of the helicity flip 
rate, and the expression has been written in terms of the Mandelstam variables. 

Given this expression we follow the “w-channel parametrization” method in Ref. 


to simplify Eq. (C4). The integration over p' is trivially performed by the 3-momentum 
conservation. The outgoing momentum k' is traded by q = k' — p, which is taken along 
the OZ axis, and a rigid rotation of dvr gives two more integrations. The momentum p is 
taken to be in the OXZ plane with an angle dpq with respect to the Z axis, and its azimuthal 
integration gives a trivial 27r. Finally, k is arbitrary in space, parametrized by its modulus 
and its polar O^q and azimuthal angle 0. 

We further introduce the energy transfer variable u, 


6{p + k — p—k')= / doj + p — k') 6 {!jJ — k + p') 


(C6) 


The two Dirac deltas in (C6) help to perform the polar integrals using 

5{u + p-k') = -dfcosOpq--^ -—^©(w-Fp), 

pq \ 2pq q J 

-k+p') = ^^s (co,e,, + ^ ^j e(i-, 

These changes bring the following conditions for the integration variables: 

p>^{q-uj) ; fc > + w) ; |a;| < g . 


Introducing all terms, we can express (|C4|) as a hve-dimensional integral: 
driR 


dt 


4 2 

e m 
2(27r)6 


-27r 


dq / doj 


dp 


dk 


'-q 


' q — io 


' q+^ 
2 


d(j) 


— q‘^ — m'^y 


F . 


(C7) 

(C8) 


(C9) 


(CIO) 


The integral in 0 can be easily performed by noting that only the s variable depends on it 
s = s(cos0). The rate is now reduced to four integrals. 


dn 


R 


4 2 

e m TT 


dt 2{27rY 


dq / du 


dp 


dk 


— g^ — Apk + 2pu — 2ku uY 




'-q 


' q — io 


' q+io 
2 


— g^ — m'^y 


q^ 


F . 

(Cll) 


The Fermi-Dirac distribution functions are expanded to linear order in the chemical 
potential so that F can be linearized in ps. 


F = 






o^P 


J' ^P{k-u) 1 g^(p+a;) _ g/3fc _ g/3p ]_ 


+ 0 


2^2 


(C12) 
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Finally, we express the solution in terms of the adimensional variables q/3 = xi ] ojjS = 
X 2 ] p(3 = X 3 ] k(3 = Xi , m(3 = y to obtain 

drifi 


dt 


1287r5 


7 ( 1 /^) 


(C13) 


where we have defined the function 


rxi 


7(r) = 


dxi 


dxo 


dx?. 


dxi 


1:1+3:2 


X 


"2 

X 2 — X^ — 4 X 3 X 4 + 2 X 3 X 2 — 2 X 4 X 2 X 2 — x\ 


qX 4 


0^3 


[x% 




y 


2\2 


qX4,-X2 _|_ 1 QX 3 +X 2 _ \ qX4, _ Y X 


Xl 


Adding a similar contribution for the left-handed fermion to get dn^/dt and using 

3 


Eq. (Cl), we obtain 


m 




(C14) 


^ 27r3^ T 

The function 7 ( 1 /^) is to be computed numerically and it acts as an IR regulator for the F 
due to the fact that it diverges in the limit ?/ —)■ 0. This function is shown in Fig. [T^ 



FIG. 13: Function 7 ( 2 /^) as a function of y = m/T. 


Fixing T = 225 MeV and a = 1/137 we hnd 


QED 



5.6 ■ 10^ fm (for m = 5 MeV) , 
6.4 ■ 10^ fm (for m = 100 MeV) , 


while for T = 500 MeV we get 


-QED _ -p-l ^ 

/ ~ ^ f - 


2.5 ■ 10^ fm (for m = 5 MeV) , 

2.6 ■ lO"^ fm (for m = 100 MeV) . 


(C15) 


(C16) 
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Evidently, the quark scattering with polarized photons cannot give any contribution to the 
helicity flip rate in the typical time scales for HICs. Results for the quark-gluon scattering) 
can be obtained from the QED result by changing a —)■ and multiplying the rate by the 
color factor 2/9 [26]. For T = 225 MeV 


QCD _ QED 




121 fm (for m = 5 MeV) , 
138 fm (for m = 100 MeV) , 


whereas for T = 500 MeV, 


(C17) 


QCD 




54 fm (for m = 5 MeV) , 

57 fm (for m = 100 MeV) . 


(CIS) 


These times have been reduced by the strong coupling constant, but they are still out of the 
limits of the lifetime of a HIC hreball, and clearly much longer than the time scales obtained 
for the change of the chirality associated to the presence chiral anomaly. 

Quark-quark scattering events might also gives a contribution to T/, but they should also 
scale as aim?/T, and although we have not explicitly computed them, we expect them to 
be also irrelevant for the change of the chiral imbalance in the HIC. 

Finally, we notice that the helicity change of a fermion (either massless or massive) by 
Landau damping is ruled out. This process can be described by cutting the fermion self¬ 
energy diagram (Fig. 14) at hnite temperature and/or density [52] (see also Ref. |53j ). 



FIG. 14: Fermion self-energy diagram, the imaginary part of which encodes Landau damping. 

The damping rate is dehned in terms of the imaginary part of the fermion self-energy 
evaluated on shell {E‘^ = + m?) 


liE) = Tr [Im S(po + ip, v){f’ + m)] 


(C19) 


Po=E 
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where the fermion self-energy can be computed in the imaginary time formalism 

d?q 


S(za;„,p) = e^T ^ j 


(2vr) 




(C20) 


with S{P—Q) and the fermion free propagator and the photon polarization function, 

respectively. 

Using the right-/left-handed helicity projectors dehned in Ref. 

1 ± 75707 ■ p 


p±(p) = 


(C21) 


one can compute the helicity-flipping damping rate as 


l{E) 


1 

4E 


Tr [Im 'P'^(p) S(po + ^7,p) V (p) {f’+ m)] 


po=E 


(C22) 


Using the algebra of Dirac matrices and the cyclic property of the trace one can easily 
check that the Dirac trace vanishes for a helicity-changing process. Therefore, Landau 
damping cannot contribute to the helicity change of a fermion. 
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